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ABSTRACT 


Context. The breakthrough developments of Cherenkov telescopes in the last decade have led to angular resolution of 0.1° and an 
unprecedented sensitivity. This has allowed the current generation of Cherenkov telescopes (H.E.S.S., MAGIC and VERITAS) to 
discover a population of supernova remnants (SNRs) radiating in very-high-energy (VHE; E > 100 GeV) y-rays. A number of those 
VHE SNRs exhibit a shell-type morphology spatially coincident with the shock front of the SNR. 

Aims. The members of this VHE shell SNR club are RX J1713.7-3946, RX J0852.0-4622, RCW 86, SN 1006, and HESS J1731-347. 
The latter two objects have been poorly studied in high-energy (HE; 0.1 < E < 100 GeV) y-rays and need to be investigated in order to 
draw the global picture of this class of SNRs and to constrain the characteristics of the underlying population of accelerated particles. 
Methods. Using 6 years of Fermi- LAT P7 reprocessed data, we studied the GeV counterpart of the SNRs HESS J1731-347 and SN 
1006. The two SNRs are not detected in the data set and given that there is no hint of detection, we do not expect any detection 
in the coming years from the SNRs. However in both cases, we derived upper limits that significantly constrain the y-ray emission 
mechanism and can rule out a standard hadronic scenario with a confidence level > 5<x. 

Results. With this Fermi analysis, we now have a complete view of the HE to VHE y-ray emission of TeV shell SNRs. All five sources 
have a hard HE photon index (T < 1.8) suggesting a common scenario where the bulk of the emission is produced by accelerated 
electrons radiating from radio to VHE y-rays through synchrotron and inverse Compton processes. In addition when correcting for 
the distance, all SNRs show a surprisingly similar y-ray luminosity supporting the idea of a common emission mechanism. While the 
y-ray emission is likely to be leptonic dominated at the scale of the whole SNR, this does not rule out efficient hadron acceleration in 
those objects. 

Key words. Astroparticle physics - Gamma-rays : general - ISM: supernova remnants - SNR : individual : HESS J1731-347, SN 
1006 


1. Introduction 


The breakthrough developments of Cherenkov telescopes in the 
last decade have enabled exploration in the very-high-energy 
(VHE; E > 100 GeV) y-ray sky with angular resolution bet¬ 
ter than 0.1° and an unprecedented sensitivity. Observations 
with the current generation of Cherenkov telescopes (H.E.S.S., 
MAGIC and VERITAS) of the Galactic plane have revealed a 
population of supernova remnants (SNRs) radiating in VHE y- 
rays. Such radiation is the signature that particles are accelerated 
to multi-TeV energies in SNRs. 


VHE observations, combined with longer wavelength infor¬ 
mation, have been used to probe the nature (hadrons or leptons) 
and the energy properties (energetics and maximum energy) of 
those accelerated particles in order to investigate the origin of 
Galactic cosmic rays (CRs). 


TeVCaflor the catalog of high-energy observations of Galactic 
SNRfQ). This sample could be described in two main categories: 
1) a population of young objects (t <5 kyrs) where the emission 
is likely associated with the SNR shell and 2) a group of older 
SNRs (t > 10 kyrs) that are mostly radiating through their inter¬ 
action with the surrounding molecular clouds (e.g. W28, IC443). 

For a number of objects in the first group with large an¬ 
gular size (Rsnr > 0.25°), the VHE emission is spatially re¬ 
solved and exhibits a shell-type morphology spatially coinci¬ 
dent with the shock front of the SNR. The members of this 
VHE shell SNR club are RX J1713.7-3946 (lAharonian et al.1 
l200l, RX J0852.0-4 622 (lAharonian etakll2007bl) . RCW 86 
(lAharonian et al.ll2009l), SN 1006 (lAcero et alJ2010l) . and HESS 
J1731-347 ( Abramowski et al1l201 ll) . In those objects the y-ray 
emission allows us to probe the population of high-energy parti¬ 
cles directly at the shock where the acceleration is taking place. 


The list of VHE sources associated with a SNR has been 
steadily growing in the last decade and a dozen detections have 
been reported so far (see, for example, the TeV online catalog 


1 http://tevcat.uchicago.edu 

2 http ://www.physics.umanitoba.ca/snr/SNRcat/, iFerrand & Safi-Harbl 

d2012h 
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Those objects also share a number of similarities. For example, 
they are evolving in low-density ambient medium (< 1 cm -3 ) and 
are bright X-ray synchrotron emitters (see Sect. I3,3l and Table [2] 
for references). 

The VHE emission can be produced either from the Inverse 
Compton (IC) scattering of electrons off the ambient photon field 
or in interactions of hadrons with ambient matter. Understanding 
the nature of the y-ray emission based solely on the VHE ob¬ 
servations is a difficult task as both mechanisms produce similar 
spectra in the VHE regime. The Large Area Telescope (LAT) on¬ 
board the Fermi Gamma-ray Space Telescope is operating in the 
high-energy (HE; 0.1 < E < 100 GeV) y-ray domain, a crucial 
energy range where the spectral signatures of the leptonic and 
hadronic scenario differ. The hadronic emission in a case with 
an F p 2 proton spectrum produces approximately an E~ 2 spec¬ 
trum in the GeV regimtQ In the case of IC emission, an E~ 2 
electron population translates into a spectrum with a slope 5 . 

Joint studies of the HE/VHE emission in a multi-wavelength 
context (including radio and X-ray observations) are providing 
new means to disentangle the different scenarios. The best ex¬ 
ample so far is the case of RX J1713.7-3946, one of the bright¬ 
est SNRs in the VHE y-ray sky and considered as one of the 
prototypes of an efficient CR accelerator. The observations with 
Fermi- LAT of RX J1713.7-3946 have revealed a hard spectrum 
at HE with a photon index r=1.50 ± 0.11 sta t llAbdo et al.ll2011 ). 
which tends to be incompatible with purely hadronic_ mod¬ 
els (e.g. lEllison et alJl2010b Izirakashvili & Ahar onianl 2010l). A 
similarly hard index (T ~ 1.4 ± 0.2 g i al ; (Yuan et al.ll2014l) was 
found for RCW 86. 

In the case of RX J0852.0-4622, HE y-rays were detected in 
spatial coincidence with the SNR and the resulting HE spectrum 
is found to be hard (although systematic errors are rather large) 
with a spectral index of T = 1.85±0.06 sta t±0.18 VVS (Tanaka et al.l 
[201 ll ). However, we note that a fraction of the HE emission seen 
by Fermi- LAT from RX J0852.0-4622 could be associated with 
the pulsar wind nebula (PWN) surrounding PSR J0855-4644, 
which would harden the spectral slope. This PWN is in spa¬ 
tial coincidence with the south-eastern shell of the SNR and 
is seen in X-rays ( Acero et~afll201 3l) and possibly VHE y-rays 
dPaz Arribas et al.ll2012). 

For HESS J1731-347 and SN 1006, no Fermi- LAT source is 
listed at the position of the two SNRs in the third Fermi catalog 
(3FGlfl lAcero et al.l20H . SN 1006 has been studied using 3.5 
years of Fermi-L AT P7V6 data bv lArava & Frutosl d2012h and no 
detection was reported. An upper limit was reported (assuming 
a point source) on a large energy range (500 MeV to 100 GeV) 
that is not well optimized to constrain the models. The last object 
of the VHE SNR club, HESS J1731-347, was studied in the HE 
domain with 3.5 years of Fermi-LAY P7V6 data by [Yang et all 
( 1201411 and no detection was reported. 

With the advent of the new Fermi-L AT reprocessed data 
(P7REP), which provide an improved sensitivity and angular 
resolution (see Sect. [2} together with a new Galactic diffuse 
model, we investigate with 6 years of data the HE counterpart 
of the two TeV shell SNRs (SN 1006 and HESS J1731-347) 
that have not yet been detected. In Sect. [2] we present the data 
analysis while in Sect. [3] we discuss the implication of the results 
with respect to the nature of the y-ray emission and investigate 
the general class properties of the TeV shell SNRs. 


2. y-ray observations with Fermi -LAT 

2.1. Data analysis 

The Fermi Large Area Telescope (LAT) is a y-ray telescope op¬ 
erating from 20 MeV to energies greater than 300 GeV. A de¬ 
scription of the instrument and of its performance is presented in 
lAtwood et al.l rt2009h . In operation since August 2008, the LAT 
provides the most sensiti ve all-sky survey in the HE y-ray regime 
(lAckermann et al.lll2012l) . 

The analysis presented here was carried out with 6 years of 
data from August 4, 2008 to August 4, 2014 using the LAT re¬ 
processed Pass 7 data (P7REP). A detailed description of the 
reprocessed data can be found in lBregeon et al.l d2013l) . This re¬ 
processing results in several improvements in the quality of LAT 
data. Among them, the point-spread function (PSF) is signifi¬ 
cantly improved above a few GeV (~ 25% smaller above 10 
GeV) and the significance of detection and precision of mea¬ 
sured photon flux is increased slightly for most sources, more 
strongly for sources with hard spectra as could be the case for 
some young SNRs. 

The P7REP SOURCE class events and the Instrument Re¬ 
sponse Functions (IRFs) P7REP_SOURCE_V15 were used for this 
study. We selected only events with energies greater than 3 
GeV as a compromise between statistics and background from 
the diffuse Galactic emission, and zenith angles smaller than 
100° to reduce contamination from the Earth limb dAbdo et al.l 
l2009f) . Time intervals when the rocking angle was more than 
52° and when the Fermi satellite was within the South At¬ 
lantic Anomaly were also excluded. The Fermi Science Too 1:0 
v9r31pl and the corresponding Galactic diffuse background 
(gll_iem_v05_revl.fit) and extragalactic isotropic background 
(iso_source_v05.txt), were used. 

In addition to the diffuse backgrounds, the model of the re¬ 
gion of interest contains the 3FGL sources within a 10° radius 
of the SNRs positions. Using this model, we computed resid¬ 
ual test statistic (TS) maps to search for additional background 
sources. The residual TS map is obtained by computing the TS 
value for an additional point source at each point of the grid in 
excess of a_given model for the region of interest (see Sect. 3.1 
of [Nolan et ali2012l . for a more detailed description of the con¬ 
cept of residual TS maps). The positions of the excesses (with a 
threshold in residual TS of 25) were used as seeds to define the 
position of additional y-ray sources not present in 3FGL. 

With 6 years of P7REP data, we investigated the GeV coun¬ 
terparts of those two SNRs using a binned likelihood fit on a 
7° x 7° region of interest. Both SNRs are approximately 0.5° in 
diameter and, for the energy range considered here, should be 
considered as spatially extended sources for Fermi-LAT. There¬ 
fore they were modeled using morphological templates derived 
from the HESS excess maps. The energy threshold for these 
excess maps is 240 GeV and 500 GeV for HESS J1731-347 
and S N 1006 respectively (lAcero et al.ll2010tlAbr amowsk i et al.l 

l2onh . 

Although SN 1006 is located at high Galactic latitude (£ = 
14.6°) and is not as affected by the Galactic diffuse emission as 
HESS J1731-347, the same energy threshold (E > 3 GeV) was 
used for both SNRs for the sake of comparison. 

2.2. Results: HESS J1731-347 


3 Assuming that the cutoff energy in the particle population is beyond 
the HE regime. 

4 The 3FGL uses 4 years of P7 reprocessed data. 


Using the model presented in the previous subsection (diffuse 
models+3FGL sources), we generated a TS map of the region to 

5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 


Article number, page 2 of[7] 



















































Acero et al.: Study of TeV shell supernova remnants at gamma-ray energies 


search for new background sources. Based on the TS map, two 
additional background point sources (PS-SE and PS-NW, see 
Fig]]} at the Galactic positions €,b = (354.34°, -1.11°), (352.93°, 
-0.30°) were detected with TS=41.19, 31.22 respectively. 

We generated the residual TS map with a model that includes 
these 2 new sources but not SNR HESS J1731-347. This resid¬ 
ual TS map is shown in Fig[T] along with the positions of the 
3FGL sources and the newly added background sources. No sig¬ 
nificant excess y-ray emission is observed in the region of the 
SNR. Within the SNR contours, the maximal TS value is 6.4. 
s To investigate the y-ray emission under the extended source 
assumption, we added the HESS template to the model using 
a power-law spectral model. The TS value of the SNR derived 
from the binned likelihood fit is 2.58/2.55/2.28 assuming a spec¬ 
tral index of 1.5, 2.0 and 2.5 respectively. 

To derive upper limits for this extended source, we used 
the HESS morphological template of the SNR only (the nearby 
source HESS J1729-345, contours shown in Fig. Q] is not in¬ 
cluded). A power-law spectral distribution with a fixed spectral 
index of 2.0 was assumed. Upper limits at 95% confidence level 
(CL) were derived using a Bayesian method implemented in the 
Python tool IntegralUpperLimit provided in the Fermi Science- 
Tools. The energy bands 3-30 GeV and 30-300 GeV were cho¬ 
sen as the best compromise between statistics and contamination 
from the diffuse Galactic emission. The resulting limits (listed in 
TableQ3 are shown in Fig.[2]together with the HESS spectral data 
points. In the most constraining energy bin (3-30 GeV), chang¬ 
ing the spectral index to 1.5 or 2.5 alters the results by 10 to 20% 
(see Table [!}. 

While the analysis was carried out for E > 3 GeV in or¬ 
der to minimize the impact of the Galactic diffuse emission, 
its contribution in the Galactic plane is non-negligible even at 
these high energies. To estimate the associated systematic un¬ 
certainty we applied the procedure used in several Fermi anal¬ 
yses of Galactic sources (e.g. the second Fermi-LAT catalog of 
y-ray pulsars, lAbdo et akl 12013b . In this procedure the best-fit 
value of the normalization of the Galactic diffuse component is 
increased/decreased by 6%0, and the flux is recomputed with the 
latter normalization parameter kept frozen. As we are dealing 
with upper limits, we decreased the normalization of the Galactic 
diffuse emission by 6% (the residuals are thus higher), froze this 
parameter, and recomputed the upper limit. The resulting upper 
limit is 1.68 x 10 12 erg cm~ 2 s~ ] in the 3-30 GeV energy range 
(a ~2% increase in comparison with the value in Table |Tj» while 
the upper limit in the 30-300 GeV energy bin was unchanged. 

2.3. Results: SN 1006 

The same procedure, described in section [2/2l is applied to the 
region of SN 1006. No background sources were detected on 
top of the standard model (diffuse emission+3FGL sources). The 
resulting residual TS map in shown in Fig. |T] and no significant 
emission is detected within the contours of the SNR. Using the 
HESS morphological template of SN 1006, the TS value from 
the likelihood fit is very small (TS=0.95). Due to the bipolar 
morphology of SN 1006, both limbs were tested independently 
by splitting the HESS template as two regions. As both regions 
are separated by -0.5°, the Fermi -LAT PSF above a few GeV is 
sufficient to disentangle the contribution from both limbs. The 


6 The specific value of 6% represents the 1.5 a deviation in the distri¬ 
bution of Galactic diffuse normalization parameters from all the fits in 
the second catalog of pulsar iAbdo et alJ2013lf . 


Index 

J1731—347 

SN1006 

SN1006 NE 

SN1006 SW 



TS values 


2.0 

2.55 

0.95 

0.13 

0.91 



Upper limits 


1.5 

1.36/6.25 

0.49/1.39 

0.19/0.55 

0.19/0.51 

2.0 

1.65/5.21 

0.47/0.97 

0.16/0.44 

0.19/0.35 

2.5 

2.05/4.32 

0.43/0.71 

0.13/0.41 

0.18/0.26 


Table 1. Fermi-LAT test statistics (TS) values obtained using the HESS 
morphological templates. For the case of SN 1006, the HESS template 
was divided into the two bright limbs to test them independently. In the 
second part of the table, the 95% confidence level upper limits (in units 
of 10- 12 erg cm -2 s -1 ) are reported for different spectral index in the two 
energy bands 3-30 GeV and 30-300 GeV. 


resulting TS from the north-eastern and south-western limbs are 
respectively 0.13 and 0.91 for a spectral index of 2.0. 

Because of the low significance of SN 1006 using 6 years 
of P7REP data, we do not expect any significant detection in the 
coming years. Upper limits using the morphological template for 
the whole SNR and for the two limbs are derived in the 3-30 
GeV and 30-300 GeV energy bands and are reported in Table Q] 
The limits for the whole SNR together with the HESS data are 
shown in Fig. [2] As SN 1006 lies 14° above the Galactic plane, 
we expect the contribution from the Galactic diffuse emission 
above 3 GeV to be negligible. We applied the same test as in 
the previous section for a systematic effect related to the Galac¬ 
tic diffuse emission normalization and the upper limits remained 
unchanged. 

3. Discussion 

To investigate the nature of the y-ray emission from both SNRs, 
we focused on the HE upper limits derived in the previous sec¬ 
tion and the VHE data collected in the literature. The data are 
then compared with a simple one-zone static model where the 
particles (electrons and protons) are described by a power-law 
with an exponential cutoff] of the form: 

dVe,p/d£ cx E~ s - p x expi-tE/El’V'f). (1) 

The sharpness of the cutoff is represented by the factor yS. In most 
of the literature, an exponential cutoff is assumed (J3- 1.0) which 
we will use throughout this study for the sake of comparison. 
However, we note that a broader cutoff with f3 = 0.5 can also 
fairly well reproduce the multi-wavelength data and in particular 
the low-energy part of the VHE spectra (see lUin et al.l[200fe 
iLi et al.f201 ll for an example on RX J1713.7-3946). 

The y-ray emission from the decay is calculated following 
the method of lHuang et al.l d2007l) where accelerated protons and 
helium nuclei collide with the interstellar medium with standard 
composition. The IC emission is derived from the scattering on 
the Cosmic Microwave Background (CMB) and the local inter¬ 
stellar infrared radiation field derived from IPorter et al . 1 (2006j). 
A lower li mit to the distance of HES S J1731 -347 of 3.2 kpc was 
derived in lAbramowski et al] d201 ll) . In the following discussion 
the distance is fixed to 3.2 kpc and therefore the electrons and 
protons energy budget (W e and W p ) should be viewed as lower 
limits. The distances to the other SNRs are listed in Table[2] 

7 We note that in all the cases, except RX J1713.7-3946, the statistics 
at HE and VHE are not sufficient to measure the shape of the cutoff. 
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Fig. 1. Residual test statistic (TS) maps for E > 3 GeV in a 4° x 4° 
region for HESS J1731-347 (top) and SN 1006 (bottom). The back¬ 
ground model is described in Sect.[2]and the SNRs are not included. The 
point sources included in the model are represented by diamonds (3FGL 
sources) and crosses (additional background sources). The source HESS 
J1729-345 seen near the SNR HESS J1731-347 is shown on the TS 
map (top) but is not included the morphological template used to derive 
upper limits. No Fermi -LAT sources are detected in the vicinity of SN 
1006. 

3.1. Constraints on emission model for HESS J1731-347 

The Fenni-LAI upper limits reported in Section [2] are shown 
in Fig. [2] with respect to the HESS data points. In a purely 
hadronic scenario both the HE and VHE emission result from 
the n° decay. In the test particle case with a proton spectrum in¬ 
dex i p =2.0 (illustrated by a red solid line in Fig. [2}, the purely 
hadronic scenario is excluded by the Fermi -LAT upper limits. In 
order to quantify to which extent we can exclude this hadronic 
scenario, we compared the log-likelihood (2i?o) obtained when 
the SNR flux is frozen to the hadronic case, to the log-likelihood 
(2z?i) given by the fit of the SNR when the flux is let free. In both 
cases the log-likelihoods are computed using Fermi data, i.e. in 



Energy (eV) 

Fig. 2. Broadband Spectral Energy Distributions (SEDs) for the SNR 
HESS J1731—347 (top) and SN 1006 (bottom). A purely leptonic (blue 
solid) and a hadronic (red solid) scenario are shown in the test particle 
case (s ep =2.0). The red dashed line shows the steepest hadronic spectral 
slope compatible with the data, respectively i p =1.5 and 1.7 for HESS 
J1731-347 and SN 1006. The HESS data poin ts for the top and bottom 
panels are extracted from Abramows ki et alJ 120 1 ll) and lAcero et all 
1201 (!) respectively and the Fermi-LAfT 95% confidence level upper lim¬ 
its are derived assuming a spectral index of 2.0. 


the 3-300 GeV range, and the photon spectral index was fixed 
to 2. For 2S?o. the SNR flux is frozen to a value of 10.21 x ICE 10 
cm -2 s~ ] in the 3-300 GeV band which corresponds to the flux 
in the purely hadronic case where s p =2.0 that connects to the 
HESS data points. 

The resulting A2z? = (££\ - .£)) is 14.68 for one additional 
degree of freedom and we therefore conclude that the hadronic 
test particle scenario is excluded at the ~5.4 cr confidence level. 
In order to accommodate the HE upper limit and the VHE data in 
a hadronic model, a proton slope s v <1.5 is required (see dashed 
line in Fig.0. 

In several aspects (age, bright X-ray synchrotron emis¬ 
sion, low-density ambient medium, VHE luminosity) HESS 
J1731 —347 is very similar to RX J1713.7-3946 (see Table 0. 
Our Fermi-LAT analysis has shown that HESS J1731-347 has 
very similar properties at HE compared to RX J1713.7-3946 
whose measured photon index is 1.5 +0.1. It is therefore likely 
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that the dominant emission mechanism in y-ray for HESS 
J1731 —347 is similar to the case of RX J1713.7-3946 (i.e. lep- 
tonic dominated). 

3.2. Constraints on emission model for SN 1006 

The upper limits derived for SN 1006 bv lArava & Frutosl(l2012[) 
were not sufficiently constraining to be able to disentangle be¬ 
tween the different emission scenarios possible. We also note 
that SN 1006 was modeled as a point source in their study. While 
this is a reasonable assumption for Fermi-LAI in the 100 MeV 
energy range, the radius of the SNR (R=0.25°) becomes non- 
negligible in comparison with the PSF above a few GeV. In con¬ 
sequence, the value of the upper limits for SN 1006 can be arti¬ 
ficially low when the point source hypothesis is used. 

Here, with 6 years of P7 reprocessed data and assuming 
an extended source hypothesis (the HESS template), we show 
in Fig. |2] the first constraining upper limits at HE. As for 
HESS J1731-347, we compared the log-likelihoods in a sce¬ 
nario where the 3-300 GeV flux is fixed to the hadronic hypoth¬ 
esis (2.19 x 10- 10 cirT 2 s -1 , red line in Fig. 0 and a scenario in 
which the flux normalization is let free. With a A2z? = 17.13, the 
hadronic scenario is excluded at a 5.8 cr confidence level. The 
steepest hadronic spectral slope compatible with the HE upper 
limit and the VHE data points gives s p < 1.7 (red dashed line in 

Fig-Q. 

When using the HESS template to model the southwest rim 
only of the SNR, we find no evidence of HE emission. While 
the synchrotron cut-off frequency is slightly higher in the north¬ 
east rim dMiceli et al.l2009l ). the southwest rim is of high interest 
to study the hadronic scenario as this region is the only one in 
the SNR where we have evidence of efficient particle accelera¬ 
tion at the shock (traced by the non-thermal X-rays), dense target 
densities (traced by the HI observations) and proof of the inter¬ 
action between the shock and the cloud (shock front is curved 
inwards in this region) as reported bv lMiceli et all (120141) . How¬ 
ever this interaction region represents a small angular fraction 
of the whole SNR and requires detailed HE and VHE spatially 
resolved spectroscopy that is at the limit of the capacities of cur¬ 
rent generation instruments. 

3.3. Source Class comparison at HE and VHE 

In order to understand the underlying emission mechanism at 
HE and VHE, we investigate and compare the y-ray properties 
of the members of the class. Among those properties, the photon 
spectral index observed at HE is a key ingredient to disentangle 
between a test particle leptonic/hadronic scenario. The HE and 
VHE spectral index of the five SNRs considered here are listed 
in Table [2] We note that in a simple one zone leptonic model, 
the slope of the electron population is in part constrained from 
the radio to X-ray synchrotron radiation. However, as there are 
usually no data between the radio and the X-ray data points, the 
slope and the synchrotron cutoff frequency can be degenerate 
parameters in particular if synchrotron cooling effects are sig¬ 
nificant and if curvature is allowed in the particle spectrum. In 
the case of RX J1713.7-3946, RX J0852.0-4622 and RCW 86 
which were detected in the Fenni-LAY data, the electron slope 
can actually be fitted over a larger energy range using GeV to 
TeV data, thus giving better leverage on the electron slope and 
reducing the impact of the aforementioned caveats. In the two 
cases where Fermi-LAI upper limits are reported, the electron 


slope can be regarded as an upper limit (i.e. the steepest slope 
allowed by the Fermi -LAT upper limits). 

For the five SNRs considered, all show, or are compatible 
with, hard HE spectral indices (1.4 < Y < 1.8) that exclude 
the standard hadronic test particle scenario. All photon indices 
(except for RX J0852.0-4622) are compatible with a test par¬ 
ticle leptonic dominated scenario where the electron slope is 
,v e = 2.0, which translates into photon spectral index of 1.5. In 
the case of RX J0852.0-4622, the slightly higher HE photon in¬ 
dex (1.85±0.06 s tat±0.18 sys t, iTanaka et al.ll20 1 ll) could be due to 
a deviation from the test particle case, a mix of hadronic and lep¬ 
tonic contributions or a possible contamination from the pulsar 
wind nebula seen around PSR J0855-4644 (lAcero et alJl2013l) 
that is located right on the south-eastern part of the SNR shell. 

Theoretical possibilities to explain a hard spectral index 
(The < 2,0) within a hadronic sce n ario include back reaction 
effects dBerezhko & Volkl 120061 IZirakash vili & Aharon ian 
20l(ih or shock-cloud interaction (llnoue et aD 120121 : 
Gabici & Aharonianll2014t) . However in young SNRs where the 
hadronic hypothesis is preferred, the measured spectral indices 
are softer or equ al to 2.0 as in Cassiopeia A (2.0 +0,l sta i + 0.1 svsf , 
Abdo et al.l2010h and Tycho (2.3±0.2 stat ±0.1 sys t, iGiordano et alJ 
20121) . 

This similarity of hard photon spectral indices in our SNR 
sample tends to point towards a common leptonic dominated 
scenario for the HE and VHE y-ray emission. We note that 
while this is probably true when looking at the spectrum aver¬ 
aged on the whole SNR, there could be some smaller subregions 
(e.g. dense clumps) where the hadronic mechanism could signif¬ 
icantly contribute to the local y-ray emission. 

Based on the model parameters compiled from the leptonic 
model fits in the literature, the SEDs in HE and VHE y-rays in 
luminosity space of the five shell SNRs are presented in Fig. [3] 
In addition to a similar HE spectral index, this comparison plot 
reveals a striking similarity in terms of peak luminosity and spec¬ 
tral shape for the SNRs considered in this sample. This similarity 
is highlighted when compared with the SED of the SNR W44 
where the evidence for hadronic emission is secure (detection 
of the 7i° decay feature, the smoking-gun evidence for hadronic 
emiss ion, in the <100 MeV energy range: Ackerman n et alJ 
120131) . The y-ray luminosity of SN 1006 in Fig.[3]is lower than 
for other SNRs. This is probably related to the SNR bipolar mor¬ 
phology and the reduced surface for efficient particle accelera¬ 
tion. If we correct for this effect by a renormalization factor of 
0.2 as discussed in iBerezhko et all (2 0091) . the peak luminosity 
is comparable to other SNRs. 

This similar y-ray luminosity in Fig. 0 was not a priori ex¬ 
pected given the fact that our sample is composed of different 
types of SN explosion and ages ranging from 1 to 6 kyrs (see 
Table |2>. Nevertheless, the sources in our sample share an im¬ 
portant characteristic that could explain part of this y-ray sim¬ 
ilarity. This characteristic is that the SNRs have evolved for 
most of their life in a low-density ambient medium (see Ta¬ 
ble 0 which allowed them to maintain a high shock velocity 
over a long period of time and therefore efficiently accelerate 
particles to high energies. In this low-density ambient medium, 
the y-ray emission is dominated by the leptonic mechanism in 
which the main source of the photon field, the CMB, is com¬ 
mon to all the SNRs. We note that this low-density argument 
is likely an over-simplification in the case of RCW 86. This 
object probably exploded in a low-density cavity blown by the 
wind of a single degenerate system, and only recently have parts 
of the shock started to interact with the border of the cavity 
(IWilliams et al.l20TTtlBroersen et al.l2014l) . As a result, the ther- 
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RX J1713 

RX J0852 

HESS J1731 

RCW 86 

SN 1006 

SN nature 

CC 

CC 

CC 

Type Ia 

Type Ia 

Distance (kpc) 

0.9-1.3 

0.6-0.9 

> 3.2 

2.3-2.8 

2.0-2.4 

Radius (pc) 

10 

12 

14 

15 

10 

Age (kyrs) 

1.6 

2-4 

2-6 

1.8 

1 

Density (cm -3 ) 

<0.02 

<0.03 

<0.01 

0.1-0.5 

<0.05 

The 

1.5 + 0.1 

1.85+0.06 

< 1.5* 

1.4 ±0.2 

<1.7* 

Fvhe 

2.32 + 0.01 

2.22 ± 0.06 

2.32 ± 0.06 

2.41 +0.16 

2.30 + 0.15 

References 

12, 17, 11,6, 1 

13,7, 18,20,9 

21,3, 22,5 

19, 10, 15, 8, 
24 

23, 2, 14, 4, 
16 

■Se 

2.15 

2.15 

2.02 

2.30 

2.10 

E cut (TeV) 

51 

25 

24 

22 

10 

W e (10 48 ergs) 

0.55 

0.38 

0.18 

1.2 

0.21 

References 

Yuan et al. 

Tanaka et al] 

Yang et al.l 

Yuan et al. 

Acero et al.l 


(2011) 

(2011) 

(201J) 

(2013) 

(2010j) 


Table 2. Summary of the physical properties of the known TeV shell morphology SNRs. The corresponding nature of the SN explosion is noted 
as core collapse (CC) or thermonuclear (Type la). The ambient medium densities shown in this table have been derived from the presence/lack 
of X-ray thermal emission. The best-fit spectral indices at HE and VHE are reported under a power-law hypothesis in the corresponding energy 
band and the errors reported are statistical. In the second part of the table, the best-fit leptonic model for each SNR is reported. The electron 
energy budget (W e ) is given for E e > 1 GeV. For RX J0852.0-4622 and SN 1006, W e was giv en for E > 10 MeV (respect ively 100 MeV) 
and the values in the table have been rescaled. Referenc es: 1) lAbdo et all <201 lh. 2) l Acero et all 120071). 3) lAcero et all d 20091). 4) lAcero et all 
(f20~ jj)A5)]Abramow slriet all d201 lh. 6) lAharonian et all <2007 alt. 7) Aharonian et all <2007bl ). 8) lAharonian et all <2009 ), 9) AllenetalJ <2015|), 
10) Bocchino_et 1- (2000), 11) ICassam-Chenai et al.1 <20041) 12) iFesen et al.) d20lt) . 13) IKatsuda et al.1 <20081+ 1 4) IKatsuda et al.1 <20091), 15 ) 


iLemoin^Gomnardetal f^2012|), 16)|MicelietalJ d 2Ql^), 17) Moriguchi et al.1 < 20051). 1 8)lPaz Arr ibas et all <201 2l). 19) ISollerman et all <20031) 
20) lTanaka et al.1 <201 lh . 21) lTian et al] <20081) . 22) lTian et all <20 1 Gd) . 23) IWinkler et all <20031) . 24) I Yuan et al.l <2014h . 


*: this work. 


mal X-rays (tracing the shocked ISM) and the non-thermal X-ray 
emission (tracing the high-energy particles) do not stem exactly 
from the same region. A similar caveat applies in the case of 
RX J1713.7-3946 a nd HESS J1731 -347, which are surrou nded 
by molecular clouds (iFukui et al.ll2012l : lFukuda et all2014h with 
possible interactions in some regions of RX J1713.7-3946. 

It is surprising that RCW 86, where part of the shock 
emits thermal X-rays tracing ambient medium densities of ~0.5 
cm -3 , shows a similar HE spectral index to RX J1713.7-3946 
where no thermal X-ray emission has been detected so far and 
where the density upper lim it is very constraining (< 0.02 cm -3 , 
ICassam-Chenai et al.ll2004l) . Because of the density difference, 
we could have expected a higher fraction of hadronic emission 
in RCW 86 and thus a steeper spectral index at HE. We empha¬ 
size that this possible hadronic contribution on small scales is 
at the sensitivity and angular resolution limit of current genera¬ 
tion instruments. We particularly anticipate the next generation 
of Cherenkov telescope CTA to carry out spatially resolved spec¬ 
troscopy on a large energy range to test those predictions. 


4. Summary and conclusion 

Using 6 years of P7REP data of the Fermi -LAT telescope, we 
have studied the HE counterparts of the SNRs HESS J1731-347 
and SN 1006. Although both objects are not detected at those 
energies, we report new upper limits that can rule out, at a 
confidence level > 5cr, a standard hadronic emission scenario 
(i p = 2.0) as the main mechanism for HE and VHE y -ray emis¬ 
sion. Given that there is no hint of detection in 6 years of data, we 
do not expect a detection with Fermi-LAT in the years to come. 

With this study, we now have a complete view of HE and 
VHE emission of the five TeV shell SNRs. All objects show 
hard spectral indices at HE (1.4 < T < 1.8) that can simply 
be explained in a standard leptonic dominated scenario. While 
the SNRs are from different type of SNe (core collapse and Type 



Energy (eV) 

Fig. 3. Intrinsic spectral energy distribution for all the members of 
the TeV shell SNR club in the GeV-TeV energy range. For the sake 
of comparison the SED from the SNR W44 (age~ 20 kyrs), where the 
evidence for hadronic emission is secure, is shown by the dashed line. 
The corresponding parameters used to produce this figure are shown 
in Table [2] except for W44 where the broken power-law model from 
lAckermann et all <20131) is used. 


Ia), have ages ranging from 1 to 6 kyrs and are evolving in differ¬ 
ent ambient media, they all show a surprisingly similar HE and 
VHE luminosity. 

We emphasize that the fact that the HE and VHE emissions 
are likely to be dominated by leptonic emission does not rule out 
efficient hadrons acceleration in those TeV shell SNRs. How¬ 
ever due to the low-density ambient medium on average, the 
hadrons do not encounter sufficiently high target densities to 
produce a level of hadronic flux that can compete with the lep- 
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tonic emission from the whole SNR. We note that in certain 
localized regions of those SNRs, the shock probably encoun¬ 
ters enhanced densities and that locally the hadronic contribution 
might become important. The next generation of Cherenkov tele¬ 
scopes such as CTA will provide the necessary angular resolu¬ 
tion to carry out detailed spatially resolved spectroscopy, which 
might unveil different emission mechanisms depending on the 
regions. We anticipate that this will bring a shift in paradigm 
from "Is the y-ray emission leptonic or hadronic dominated ?" to 
"In which region of the SNR is the emission leptonic/hadronic 
dominated ?". This evolution might be comparable to the ques¬ 
tion about the nature of the X-ray emission of SN 1006 in the 
80’s (see iBecker et al.lll98(i and reference therein) which de¬ 
bated whether the emission was Crab-like (power-law) or Tycho- 
like (thermal). The ASCA satellite with an improved energy 
range, sensitivity, and angular resolution revealed that both emis¬ 
sion mechanisms co-exist but with a different spatial distribution 
dKovama et al.lll995h . 
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